Spin correlations in the two-leg antiferromagnetic ladder in a magnetic field 
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We study the ground-state spin correlations in the gap less incommensurate regime of a S — 1/2 
XX Z chain and a two-leg antiferromagnetic ladder under a magnetic field, in which the gapless ex- 
citations form a Tomonaga-Luttinger (TL) liquid. We calculate numerically the two-spin correlation 
functions and the local magnetization in the two models using the density-matrix renormalization- 
group method. By fitting the numerical results for an open XX Z chain of 100 spins to correlation 
functions of a Gaussian model, we determine the TL-liquid parameter K and the amplitudes of the 
correlation functions. The value of K estimated from the fits is in excellent agreement with the exact 
value obtained from the Bethe ansatz. We apply the same method to the open ladder consisting of 
200 spins and determine the dependence of K on the magnetization M. The K-M relation changes 
drastically depending on the ratio of the coupling constants in the leg and rung directions. We also 
discuss implications of these results to experiments on the nuclear spin relaxation rate 1/Ti and 
dynamical spin structure factors. 
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I. INTRODUCTION 

Spin ladder systems have been studied extensively over 
the past decade U There are reasons why the ladders 
have attracted so much attention. Firstly, they natu- 
rally interpolate one- and two-dimensional systems and 
may provide some hints to better understand the high- 
temperature superconductivity which occurs in square 
lattice CuC>2 planes. Secondly, spin ladders themselves 
have interesting physics and deserve through investiga- 
tion in their own right. One of their most surprising 
properties is that low-energy physics of spin ladders de- 
pends drastically on the number of legs. Spin- 1/2 antifer- 
romagnetic (AF) ladders, for example, have a finite gap in 
the spin excitation spectrum in the even-leg case, whereas 
they have no gap in the odd-leg case. The ground state of 
an even-leg ladder is a spin singlet and its properties can 
be understood from the short-range resonating-valence- 
bond picture.B This spin-gap behavior has been observed 
experimentally on S = 1/2 two-leg ladder compounds pa 
such as SrCu203 and Cu2(C5Hi 2 N 2 )2Cl4. 

A gapless phase can appear in even-leg ladders when 
an external field h is applied. If the field h is larger than 
a critical field h&, which is equal to the spin gap, and 
if it is smaller than the saturation field h C 2, then the 
ground state has a nonzero magnetization M and the en- 
ergy gap between the ground state and the first-excited 
states vanishes. The gapless mode has been shown to be 
described as a Tomonaga-Luttinger ffJj) liquid both in 
the strong- and weak-coupling limit ,El~Q where the cou- 
pling in the rung-direction J± is much larger or much 
smaller than the one in the leg-direction J|| , respectively. 
However, the M dependence of the TL-liquid parame- 



ter K , which governs the spin correlations in long wave 
length, has been obtained analytically only in the strong- 
coupling limit and it remains as a nontrivial problem to 
determine K for general J±/J\\ . In the gapless phase the 
system shows incommensurate spin correlations since the 
Fermi wavenumber of Jordan- Wigner fermions is shifted 
from 7r/2 in the presence of the magnetic field which acts 
as a chemical potential for the fermions. The wavenum- 
ber Q characterizing the incommensurability of the gap- 
less mode varies continuously as h increases. This gapless 
incommensurate (IC) phase is in fact in the same univer- 
sality class as the one-dimensional S = 1/2 XX Z model 
in a magnetic field, as we will see. 

In this paper, we study low-energy properties of the 
S = 1/2 two- leg AF ladder in the gapless IC regime 
for broad range of Jj_/Jn. We show that the system 
is a TL liquid for arbitrary J±/J\\ and determine the 
M dependence of K numerically. To this end, we com- 
pute numerically the ground-state spin-correlation func- 
tions and the local magnetization in the open ladders us- 
ing the.-density-matrix renormalization-group (DMRG) 
methodclll3 and extract the TL-liquid parameter by fit- 
ting the data to correlation functions obtained from the 
Abelian bosoijdzation. This method was applied in our 
previous workl to the S — 1/2 XX Z chain at h = and 
proved to be effective in determining both the TL-liquid 
parameter and amplitudes of correlation functions. In 
order to demonstrate the validity of the analysis in the 
gapless IC phase, we first apply it to the S — 1/2 XX Z 
chain for h > 0. The model is exactly solvable by the 
Bethe ansatz and the TL-liquid parameter K can be cal- 
culated for arbitrary value of M. It thus provides a good 
test ground to check accuracy of our method. We find 
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that K estimated from the DMRG data is in excellent 
agreement with the exact calculation. We then apply 
the same method to the two-leg ladders in a magnetic 
field. Our numerical data of correlation functions are fit- 
ted well for broad range of J±/J\\ to the formulas based 
on the bosonization approach, confirming that the gap- 
less modes are in fact in the universality class of a TL 
liquid for arbitrary J±/J\\. The M dependence of K ob- 
tained in the large J±/J\\ limit agrees with the analytic 
result obtained through mapping to the XXZ chain. In 
this limit K is less than 1 for < M < 1. As J±/J\\ 
decreases, K increases and become larger than 1 for in- 
termediate values of M. Our numerical result indicates 
that K takes a universal value 1 for any J±/J\\ in the 
limits M -> and M — ► 1. 

The plan of the paper is as follows. We first review 
the Abelian bosonization approach to the S = 1/2 XXZ 
chain under a magnetic field in Sec. II A. The formulas 
of the spin correlations and the local magnetization in 
finite open chains are presented. In Sec. II B, we show 
numerical data for the XXZ chain of L = 100 sites ob- 
tained from the DMRG calculation and fit the data to the 
functions given in Sec. II A. In Sec. Ill A, we briefly re- 
view some relevant results of the previous analytic studies 
on the two-leg ladders in the strong- and weak-coupling 
limits. The DMRG data and the results of fitting on the 
open ladders with L — 200 sites are shown in Sec. Ill 
B. The M dependence of K for various values of J±/J\\ 
is obtained. Its implications to NMR and neutron scat- 
tering experiments are briefly mentioned. Finally, our 
results are summarized in Sec. IV. 



II. XXZ CHAIN 

A. Bosonization approach 

In this section, we consider spin- 1/2 XXZ chains with 
open ends in a magnetic field h. The Hamiltonian is 
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ch 



7~Lo + 7~Lh 
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with 



L-l 



1=1 
L 



1=1 



where Si are S — 1/2 spin operators and (Si,S 



l' A 



SfSf, + SfSf, + ASfSf,. We assume the system size L 
to be even throughout this paper and treat only the case 
where J > and < A < 1. We note that the Hamilto- 
nian (|l|) for — 1 < A < 1 can bej-splved exactly by Bcthc 
ansatz for arbitrary values of foJilnij 



We use the standard Abelian bosonization techniques 
to analyze spin-spin correlation functions at zero tem- 
perature. We basically follow the scheme presented in 
Rcf. |l4| and generalize it to the case of open chains in 
magnetic fields. The bosonization formulas in the ab- 
sence of magnetic fields are reported in Ref. |§|. 

The low-energy dynamics of XXZ chains is described 
by the Gaussian model, Ej 
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ch 



L+l 



dx 



dx 



dx 



(2) 



where v is the spin-wave velocity. The continuous vari- 
able x is identified with the site index I under the as- 
sumption that the lattice spacing equals unity. The 
bosonic fields <p{x) and 4>{x) obey the commutation re- 
lation [4>{x), (j>(y)] = —i@(x — y), where <d(x) is the step 
function. The spin operators in the original Hamiltonian 
(0) are related to the bosonic fields by the relations, 



, - J-^ + «(-!)< sin 
1 2nR dx y ' V R J 



Sf = 



Sf = exp[-i2irR(f>(l)} 



^'M^ ) + c (~ 1 ) i 



(3) 
(4) 



with a, b, and c being real constants. The param- 
eter R determines the exponents of correlation func- 
tions. We also introduce the TL-liquid parameter K by 
K = 1/{AttR 2 ); K = 1 in the XY case (A = 0), and 
K = 1/2 in the Heisenberg case (A = 1) at h = 0. From 
Eq. (||), Sf is written as 



Sf 



c(-l) i cos[2 7 ri?0(O] 



ibsm[2%R4>(l)]sm (Jjp) ■ 

(5) 



The second term with the coefficient ib is Hermitian due 
to the commutation relation [(f)(1), <f>(l)] = —i/2. The 
open boundary conditions are translated to the boundary 
conditions on the bpsonic fields at the two phantom sites 
I = and I = L+lM(f>(0) = and <p(L+l) = 2TcRLM ch , 
where M c h is the magnetization per site, 



1 L 
M ch = -J2Sf 



(6) 



1=1 



The total magnetization LM c h is an integer for even L. 
These boundary conditions lead to the mode expansion, 

Hx) = -^00 + £ (a n + 4) , (7) 

71—1 

(an-alr), (8) 



where q n = 7rn/(i+l), [ipo, <po] = i, and boson an- 

nihilation operators obeying ^,0],] = 5 mn . Note that 



2 



the commutation relation between <p(x) and 4>{y) men- 
tioned above is satisfied. The lowest energy state |M c h) 
in the subspace in which the magnetization per spin is 
M c h is a vacuum of a n 

a n \M ch )=0 (9) 

and an eigenstate of (f>o 

(f>o\M ch ) = 2nRLM ch \M ch ) . (10) 

We may regard 0o as a coordinate variable along a ficti- 
tious ring of radius 1/2ttR and take <po = —id/d<f>o to be 

(S?S?,) 



(S?) 



its momentum conjugate. The state |M c h) is then pro- 
portional to exp(i2TrRLM c ii<f ) a)- The bosonization for- 
mulas (||) and (U) represent only the leading contribu- 
tions. In the next order S[~ has a term of the form 
(-l) l exp[-2niR4>(l)]cos[2<f>(l)/R]. We will, however, ig- 
nore this contribution because it yields only subleading 
corrections that disappear quickly for large |Z — Z'|. 

Using Eqs. one can evaluate the two-spin cor- 

relation functions {SfSp) (a = x, z) and the local mag- 
netization (Sf ) in open chains, where (• • •) denotes the 
expectation value in the state |M c h). Brief account of 
their derivation is given in Appendix. Here we present 
the final results: 
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The wavenumber q, characterizing the IC character of the 
spin correlations in a magnetic field, is related to M c h by 



2-!rM ch L 
L + l 



(17) 



The factor L/(L + 1) appears as a result of the open 
boundary conditions. Under the periodic boundary con- 
ditions q should be simply equal to 2nMch, because the 
first term in Eq. (ra) is 4>qx/L in this case. This term 
must be 4>qx/(L + 1) in the open-boundary case in order 
for <j>(x) and 4>(x) to satisfy the commutation relation in 
the interval [0, L + l]. We emphasize that Eqs. ( |l2| ) and 
( |l3| ) reproduce the exact results for the XY chain when 



rj = 1/2 and a = — 1/n. As is well known, the XX Z 
spin chain is equivalent to a model of spinless fermions 
with nearest-neighbor interaction. In this model, Sf is 
none but the fermion density, and the oscillating term in 
Eq. (|l^) corresponds to the Friedel oscillations near the 
open endsllalla 

In the thermodynamic limit (L — > oo) with \l — L/2\ -C 
L and |Z' — L/2| <C L, the spin correlations have the 
asymptotic forms 



(SfS?,) 

(Sfsp) 



= A. 



(-1) 



i-v 



= M ch 2 + A z 
1 



A 



cos [q(l - I')] 
: |Z - l>\v+Vv ' 

, .si-j> cos [q(l-l')] 
[ 1 jZ-Z'l 1 /') 



4tt 2 77|Z -Z'| 



(18) 



(19) 



where the correlation amplitudes A x , A x , and A z are re- 



lated to the numerical constants a, 6, and c by A x 



V2, 



-Arc = b 2 /4, and A 2 = a 2 /2. We can therefore estimate 
the TL-liquid parameter and the correlation amplitudes 
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in the thermodynamic limit by extracting the fitting pa- 
rameters i?, a, b, and c from the numerical data on a 
finite system with use of Eqs. @, (|l2|), and |f3|). At 
the same time, the TL-liquid parameter K can be calcu- 
lated exactly for any M c h by solvipg_an integral equation 
obtained from the Bethe ansatzE3'E!l We will compare 
our estimates of K obtained from the fitting procedure 
with the Bethe ansatz results in the next subsection. 



B. Numerical results 



Using the DMRG method,BlI3 we computed the two- 
spin correlation functions (SfSf;) (a = x, z) and the 
local magnetization (Sf) in the L = 100 open chains. 
The two-point functions were calculated for I = ro — r/2 
and I' = ro + r/2, where ro = L/2 for even r and 
r = (L + 1) /2 for odd r. The calculation was performed 
for each lowest-energy state of Ho in the subspace of var- 
ious values of M c h- We employed the finite system al- 
gorithm of improved version.Ej The maximum number of 
kept states m is 100. We estimate the numerical error due 
to the truncation of the Hilbert space from the difference 
between the data with m = 100 and those with m = 70. 
The estimated errors for (SfSf,), (SfSf), and (Sf) are, 
at largest, of order 10 -5 ,10 -6 , and 10 -6 , respectively. 

In Fig. [I], we show the spin correlations (S"Sj?) (a = 
x, z) and the local magnetization (Sf ) at A = 0.5 for 
three different values of M c h- The DMRG data are shown 
by open symbols whose sizes are larger than the trun- 
cation error mentioned above. Taking R, a, b, and c as 
fitting parameters, we fit the numerical data to Eqs. (|l"l| ) 
(H). The results of the fitting using the DMRG data of 



(SfSft) for 10 < r < 90 and of (Sf) for 1 < I < 100 
are also plotted in the figure by the small solid symbols. 
One can see that the fits are in excellent agreement with 
the DMRG data, proving the validity of the bosonization 
formulas 

For various values of M c h and A, we determined the pa- 
rameters R, a, b, and c by the fitting procedure. In doing 
so we used numerical data of several ranges, 10 < r < 80, 
10 < r < 90, 20 < r < 80, and 20 < r < 90 for 
the two-spin correlation functions and 1 < I < 100 and 
10 < I < 90 for the local magnetization. We take the 
mean and the variance of the fitting parameters obtained 
for the different ranges of r and I as the estimated value 
and the error of the estimates, respectively. The TL- 
liquid parameter K = 1/(AttR 2 ) estimated from the nu- 
merical data of (SfSf,) is plotted as a function of M c h in 
Fig. ||. The, exact values obtained from the Bethe ansatz 
methodE3c3 are also shown as dotted lines. The agree- 
ment is excellent. We also estimated K from the fitting 
of (SfSf) and (Sf ) and obtained similar results as Fig. ||. 
We found, however, that the estimates from the last two 
correlators show some deviations from the Bethe ansatz 
results when K is small. We do not exactly know why 
they deviate. One possible reason might be the effect 
of the leading irrelevant operator neglected in the Gaus- 
sian model that becomes marginal at K = 1/2. In the 
XY regime of our interest, spins have stronger correla- 
tions in the S x and S y components than in S z , and thus 
we may expect that (SfSf) should give us most reliable 
estimates. 

For the correlation amplitudes A x and A z , Lukyanov 
and Zamolodchiko v c on jectured the exact formulas which 
are valid at ft, = 0,E3c3 



.4 



LZ 



1 



r( 



2(1— tj) ^ 



V^(2(TH)) 
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sinh(r/£) 



A L - — 



TV ^ \ 
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1/7, 



exp 



dt 

T 



sinh(i) cosh[(l — rf)t] 
sinh[(2?7 - l)t] 2-q 



r/e 



-21 



sinh(?7i) cosh[(l — rj)t] 



<1 



(20) 
(21) 



where T(x) is the Gamma function. Equations j20| ) 
and (|2l]) have been confirmed numerically.Sa In Ta- 
ble @, we give our estimates of the correlation amplitudes 
A x = c 2 /2 and A z = a 2 /2 obtained from the fitting 
of (SfSf,) and (Sfl for < M ch < 0.5, together with 
the exact values ( po|) and ( pl| ) at Af c h = 0E3 As can be 
seen in Table | (a) , A x decreases monotonically from the 
value given by Eq. (|2(]) to zero as M ch increases from 
to 1/2. Thus, A x depends not only on K but also 
on M c h. (See, for example, the data for A = where 
K takes a constant value 1 for any M c h.) When M c h 
approaches 1/2, where K — > 1, A x seems to go to zero 
linearly for any A. This can be easily understood once 
we consider one-magnon contribution to the correlation 
function. On the other hand, A z decreases monoton- 
ically from the number given by Eq. (^lj) to the uni- 



versal value A~-.= 1/(2tt 2 ) ~ 0.05066 as M c h increases 
from to l/2.t3 An exception is the case A = 0, where 
A z = l/(27r 2 ) for any M c h- The convergence of A z to 
l/(27r 2 ) at M c h — * 1/2 is consistent with the fact that 
the correlator (SfSf) must take a constant value 1/4 at 
M ch = 1/2. The right-hand side of Eq. @ equals M ch 2 
when 77 = 1/2, q = 7r, and A z = 1/(2tt 2 ). 

III. TWO-LEG AF LADDER 

Encouraged by the success in the last section, we study 
the two-leg AF ladders in a magnetic field using the same 
method. We begin with a brief review of the analytic 
results on the ladder in the strong- and weak-coupling 
limits. 
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A. Review of Analytic Results 

The Hamiltonian of the open two-leg ladder studied in 
this section is given by 

L-l L 

H = J\\ 2_j / ,QS>,i) S^i + i)a + Jx / ](Si,i, S2,i)a 
11=1,2 ;=i 1=1 

11=1,2 1=1 

The anisotropy A is introduced for generality. We as- 
sume that the coupling in the leg- and rung-direction, J|| 
and Jj_, are positive (antiferromagnetic). The spin ladder 
has an excitation gap in weak magnetic fields h < h c \. 
We concentrate on the ladder in the gapless IC regime, 
i.e., in the case h c i < h < h c %. We denote the ratio 
J±/J\\ by j hereafter. 

We begin with the strong-coupling limit (j 3> 1), for 
which a simple intuitive picture is available. It is known 
that the system in this limit. rap. be mapped to an ef- 
fective 5 = 1/2 XX Z chainJO'ESEZI as we explain below. 
Let us first assume J\\ = 0. In this case, an eigenstate of 
TL is written as a direct product of rung states. At each 
rung two spins S\j and S^,; are either in a singlet state 
\si) = (| TJ-) — I 4t))/\/2 or in one of the triplet states, 
> = I TT), |t?) = (I U) + I 4T))/V2, and \tj) = | [[). 
When h is small, the ground state consists of a product 
of the singlet rungs. As the field h increases, the energy 
of the state \tf) becomes lower, and at h = Jj_(l + A)/2, 
the state degenerates with \si). We can thus analyze the 
low-energy properties of the system for h~ Jj_(l + A)/2 
by retaining only the two lowest-energy states \si) and 
|f^~) for each rung. We may regard the two states as 
"down" and "up" states of an effective S = 1/2 spin, 



5Z QZ 



J 2,l 



$2,1)1 
1 

2' 



Sf 



1 

V2 V 



Sir) 



(23) 
(24) 



The effective spins Si are the only low-energy degrees of 
freedom, and their dynamics is governed by 



L-l 



A 



H — J|| ^2 (&h Si+i)a/2 ~ ~rJ\\ ( <$! + S 



1 = 1 



1 + A A \ x-^ ~ 

h — —Jx - — J\\ ) 22 s i + const - ( 25 ) 

1=1 



Note that the anisotropy of the effective 5 = 1/2 XX Z 
chain is a half of the anisotropy of the original lad- 
der, A/ 2. Besides the bulk effective magnetic field 
h — (1 + A)Jj_/2 — AJii/2, there is an additional field 

— AJ||/4 applied only to the boundary spins Sf and 5£. 
It induces oscillating magnetization near the boundaries 
superposed on the Friedel oscillation which is already 
present at any M 7^ without the boundary field. The 
effect of the boundary field may be cancelled by adding 
an extra term 



H' = ti 



E 

11=1,2 



11,1 



s 



(27) 



to the original ladder Hamiltonian (|22j), where hi 
AJ||/4 for j ^> 1. Now we define 



Sq,i — Si,l + Sz,i, 
S n ,l = Si,i — $2,1- 



(28) 
(29) 



From the mapping explained above, we conclude that the 
two-spin correlations (S% t S% t ,} and (Sq ; Sq t ,) and the lo- 
cal magnetization (Sq ; ) in the open ladder H + Tt' in the 
limit j > 1 are given by the corresponding correlators in 
the XX Z chain. We thus obtain 



(Sl l Sl l ,)=2X(l,l';Q) 

(SqiSqii) = 



(30) 



1 ~z\ \ ~z 



\ + \[z(l;Q) + z(l';Q)] + Z(l,l';Q), (31) 



(Sh) = l - + z(l,Q), 
where the wavenumber is 



27ri fnr 1 

q= ltt r-2 



(32) 



(33) 



In the limit L — ► 00 the two-spin correlation functions 
reduce to 

f_l)l-l' 

/qx qx \ _ O A V I 

ZAx( - L > |J _ y\2K+{l/2K) ' ^ 

WW = M2 - 4n2 v \l-l'\* 

cos[2itM(1 - I')} 



+A Z 



l\2K 



(35) 



This mapping is derived in lowest order in Jn/Jx an d 
valid for the entire IC region of < M < 1, where M is 
magnetization per rung 



1 L 



(26) 



1=1 



Note that they can be obtained from Eqs. (y_8|) and (|19| ) 
by replacing q and M c h with 2tt(M— 1/2) and M, respec- 
tively. On the other hand, the correlations (Sq ;SJ;/) and 
(S* iS* v ) decay exponentially because Sq ; and S^ ; al- 
ways create the high-energy rung states \tf) and \t7) as 
a virtual excited state. 
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Next, we consider the opposite case, the weak-coupling 
limit (j <C 1). The system in this limit has be^n-, in- 
vestigated with the Abelian bosonization method.ErQ In 
these studies, two chains are first bosonized indepen- 
dently, and then the interchain coupling J± is treated 
perturbatively.Ej Four bosonic fields <f>± (x) and 4>± (x) are 



introduced, where 



and 



[<j>— and <j>J) are the sym- 
metric (antisymmetric) combinations ot,bosonic fields of 
each chain. All the fields are massivea when h < h c \. 
In the IC regime of h c \ < h < h C 2t on the other hand, 
the fields 0+ and </>+ become massless while the fields <f>- 
and 4>- remain massive. The low-energy effective Hamil- 
tonian for the gapless modes has the same form as that 
of the S =1/2 XX Z chain, Eq. (§). Furthermore, the 
spin correlation functions (S% t S% y) and (Sq;Sq;,) have 
the same r dependence as in the strong-coi»ling limit 
(but with different values of K, a, b, and c).LrE3 The coiu 
relators (Sq;Sq;/) and (S* t S% t ,) decay exponentially,LJ 
because they involve the massive fields. This result also 
matches the strong-coupling limit. Moreover, the incom- 
mensurate wavenumber for the short-ranged correlators 
is q = nM , which is different from the IC wavenumber 
for the quasi-long-ranged correlators q = 2ir(M — 1/2). 
For example, it was found thatQ 

(S" S z )-A( -U-ri/£ Cos[7rM(/ - 1')] 

V^r^TT.W - A *\ l ) e ij _ ^11/2+1/477 ' 



where £ is a correlation length for a massive mode and 
A z is a constant. 

We have seen that, both in the strong- and weak- 
coupling limits, the low-energy physics of the two-leg lad- 
der in the gapless regime is in the same universality class 
as the XX Z chain in a magnetic field. In particular, the 
spin correlation functions (S^,S^;/) and (Sq;Sq Z ,) and 
the local magnetization (Sq t ) have the same forms as the 
corresponding functions in the XX Z chain, but with the 
shifted wavenumber r 2jr(Af — 1/2) and with different val- 
ues of K, a, &, and c.E3 It is then very natural to postulate 
that the universality is not restricted to the two limits but 
holds for any j. This allows us to use Eqs. (|30|)-(|32|) for 
analyzing the correlation functions in the ladder for any 
j and M (0 < M < 1). We can thus determine the TL- 
liquid parameter K of the ladder by fitting the numerical 
data to Eqs. (^)-(^) in the same way as we did for the 
XXZ chain. The result is presented in the next subsec- 
tion. Finally, we may expect that the \l — 1'\ dependence 
of the short-ranged correlators (SqjSq^) and (S* t S* t ,) 
obtained in the weak-coupling analysis, such as Eq. (36), 
should also be valid for any j and M (0 < M < 1). 



B. Numerical Results 

Here we present the result of the DMRG calculation of 
the two-spin correlations (S^ ; S^;,) and (Sq;Sq;,) and 
the local magnetization (Sq;) in the ladder. Taking R, 



a, b, and c as free parameters, we fit the data to the for- 
mulas (|30l)-(32) for several values of j and estimate the 
M dependence of K = 1/(AttR 2 ). The numerical cal- 
culations were performed for the open ladder 7i + 7i' of 
L = 100 rungs (200 sites) for j = 10.0,2.0,1.0,0.5 and 
A = 1.0,0.5,0.0 using the finite-system DMRG method 
of the improved version. We calculated the two-spin cor- 
relations for / = ro— r/2 and I' = ?'o+r/2, where tq = L/2 
for even r and rg = (L + l)/2 for odd r. The maximum 
value of the kept states m is 160. From the difference 
between the data with m — 160 and those with m = 120, 
we estimate the numerical error due to the truncation. 
The estimated errors for (S^S^,), (SfoSfo) , and (Sf) 
are, at largest, of order 10~ 4 , 10 -6 , and 10 -6 , which is 
almost negligible. In the course of the calculation, we 
optimized the value of the extra boundary field h' to 
minimize the effect of the boundary field.EiJ For finite j, 
however, the boundary effect cannot be eliminated com- 
pletely since it can be represented by the form of TC only 
in the strong-coupling limit j> 1. As a result, the two- 
spin correlation functions and the local magnetization in 
the open ladder H + H' might deviate from the expected 
form, Eqs. (^0|)-(p2[), near the boundaries. For this rea- 
son we used data of smaller range of r for the fitting than 
in Sec. II to reduce the unwanted boundary effect. We 
chose the regions 10 < r < 70, 10 < r < 80, 20 < r < 70, 
and 20 < r < 80 for the fitting of the two-spin corre- 
lators and 10 < I < 90 and 20 < I < 80 for the local 
magnetization. As in Sec. II, we regard the mean and 
the variance of the fitting parameters obtained for these 
different ranges of r and I as the estimated value and 
the error of the estimates, respectively. Incidentally, we 
have also checked for (j,M) = (10.0,0.5) that, without 
the boundary field h' , the local magnetization (Sq;) has 
the Friedel oscillations induced by the effective boundary 
field. We found that the oscillations decay algebraically 
into the bulk with thej-exponent K, as expected from the 
bosonization analysis .E3 

The numerical data of (S* ;S* ;,), (S 2 ; S 2 ;,), and (S 2 ,) 
for j = 10.0 and 1.0 with A = 1.0 (Heisenberg case) are 
shown in Figs. [3] and ^ by open symbols, whose sizes 
are larger than the truncation error mentioned above. 
The small solid symbols in the figures are the fits to the 
DMRG data of the two-spin correlations for 10 < r < 80 
and to those of the local magnetization for 10 < I < 90, 
respectively. It is clearly seen that the fitting works ex- 
tremely well for j = 10.0, confirming the validity of the 
formulas (|30|), j3l|), and (32). Furthermore, the agree- 
ment between the numerical data and the fits at j = 1.0 
is also quite good except some deviations near the bound- 
ary, indicating that the formulas are accurate for the 
intermediate-coupling regime of j as well. We note that 
the quality of the fitting is also good for other values 
of j and A that we have examined. We therefore con- 
clude that the gapless mode of the two-leg ladders is a 
TL liquid for arbitrary j, and accordingly, the properties 
of the strong- and weak-coupling ladders are smoothly 
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connected. 

Next we show in Fig. || the M dependence of K esti- 
mated from the data of (S* v ) for various j in both 
the HeisenbeM (A = 1.0) and XY (A = 0) cases. In the 
earlier studyE3 the exponent r\ was obtained for j = 5.0 
in the Heisenberg case only. We note that the esti- 
mation from (S% t S% is more reliable than that from 
(SqjiSqjz) oy J.Sqi) as we have seen in the XXZ chain. 
Theoreticallju it is expected that K should approach the 
universal value K = 1 when M — > as well as when 
M — > 1, since the system is equivalent to the dilute limit 
of hard-core bosons. Although the data for M — > have 
large error bars, we may conclude that our results for 
the Heisenberg case are consistent with the theoretical 
prediction. Our results for the XY case show a more 
subtle feature. At first sight the results for weaker cou- 
plings (j — 1.0 and 0.5) do not seem to approach K = 1 
as M — > 0. We think, however, that K changes very 
rapidly at small M to approach K = 1, in view of the 
data for j — 10.0 and 2.0, which are consistent with the 
theory. Unfortunately, it is difficult to numerically esti- 
mate K for small M with high accuracy to resolve this 
issue. 

In the strong-coupling limit j ^> 1, the ladder system 
with anisotropy A is equivalent to the 5=1/2 XX Z 
chain with anisotropy A/2, as explained in the previous 
subsection. Figure [| clearly shows that the estimated 
value of K for j = 10.0 is consistent with the anticipated 
behavior shown as the dotted curves in both the Heisen- 
berg and XY cases. As j decreases, K increases mono- 
tonically for any M (0 < M < 1). Thus, K is always 
larger than 1 in the XY ladder because K — > 1 for any 
M in the large j limit. In the Heisenberg ladder, on the 
other hand, K is smaller than 1 in the strong-coupling 
limit, as expected from the mapping to the XXZ chain. 
Upon decreasing the interchain coupling j, K starts to 
increase and the K-M relation changes from a concave 
curve to a convex one. We note that the similar behavior 
is observed also in the ladder with A = 0.5: As j de- 
creases, the K-M relation changes from a concave curve 
at j 3> 1, corresponding to the behavior of the XXZ 
chain with anisotropy A/2 = 0.25, to a convex one. We 
thus consider that this behavior of the K-M curve is an 
universal feature for < A < 1. 

The TL-liquid parameter K determines the long- 
distance behavior of correlation functions in the thermo- 
dynamic limit. For example, the leading term of the cor- 
relator (S§ tl Sfo,)-M 2 decays as cos[2nM(l-l')]/\l~l'\ 2K 
for K < I, while it decays like \l - l'\~ 2 for K > 1. 
Hence, in the ladder with < A < 1 the leading term of 
the correlator changes from cos[2irM(l — — l'\ 2K to 
\l — l'\~ 2 at a critical value j c {M) as j decreases, while 
in the XY ladder the leading term is always \l — l'\~ 2 . 
On the other hand, the correlator (S% v ) decays as 

(— l)' - ' /\l — l'\ 1 / 2K in the whole range of K covered 
in Figs. ^| (a) and || (b). The temperature dependence 
of the spin-lattice relaxation rate 1/T\ in NMR experi- 



ments is directly related to the TL-liquid parameter K 
through the decay exponent of the most slowly decay- 
ing correlation.!! From the behavior of K-M relation ob- 
tained above, we find that the correlator (S% v ) de- 
cays most slowly for any j, M, and < A < 1. We there- 
fore conclude that at low temperatures the relaxation 
rate of the ladder in the gapless regime always shows a 
power-law divergence l/Tj oc T~ 1+ ( X / 2K \ 

Figure H shows the numerical result of (S* ,S* v ) for 
the Heisenberg ladder at j = 0.5. It exhibits exponen- 
tially decaying oscillatory behavior. From the period A of 
oscillations, we obtain the IC wavenumber q — 2ir/\ as a 
function of M; see the inset figure. The result confirms 
the theoretical prediction q = nM. This IC wavenum- 
ber q tells us that the massfwe magnon dispersion has a 
minimum excitation energy atl3 q = n—q = 7r(l— M). Ac- 
cordingly, the dynamical spin structure factor S% z (q, u) 
should have a power-law divergence along the energy dis- 
persion which is roughly shifted by irM from that of the 
triplet magnon dispersion in the absence of the magnetic 
field. It would be interesting if this feature is observed 
by inelastic neutron scattering experiments. 



IV. CONCLUSIONS 

In this paper we have studied the ground-state spin 
correlations in the gapless IC regime of the S = 1/2 
XXZ chain and the two-leg AF ladder in a magnetic 
field. We have used the 5=1/2 XXZ chain as a first 
test ground to apply the method we developed in our 
previous work: We numerically computed the two-spin 
correlation functions and the local magnetization by the 
DMRG method and fit the results to functions which are 
obtained using the bosonization technique. The fitting 
parameters are the TL-liquid parameter K and the am- 
plitudes of bosonic operators. We found good agreement 
between K estimated from the fitting and K calculated 
from the Bethe ansatz. As a byproduct we obtained the 
amplitudes of the dominant terms in (SfSf,) and (SfSf,). 

We have applied the same technique to the two-leg AF 
ladder in the gapless IC regime. It has been known that 
in both the strong- and weak-coupling limits the low- 
energy excitations in the ladder are regarded as a TL 
liquid like the XXZ chain in a field. We fit our DMRG 
data of the two-spin correlation functions and the local 
magnetization of the ladder to the same bosonization for- 
mulas we used in the analysis of the XXZ chain. The 
fitting worked very well not only in the strong- and weak- 
coupling limits but for broad range of the interchain cou- 
pling strength j. We thereby confirmed that the low- 
energy gapless excitations are indeed described as the 
TL liquid for any j and the properties of the strong- and 
weak-coupling ladders are smoothly connected. For sev- 
eral values of j, we have determined K, which shows non- 
trivial j and M dependences (Fig. ||). It turned out that, 
for any M(0<-M<l),irisa monotonically decreasing 
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function of j . In the ladder with anisotropy < A < 1 
the K-M relation changes from a concave curve at j 1 
to a convex one as j decreases, while in the XY ladder 
(A = 0) it changes from a line K — 1 at j ^> 1 to a con- 
vex curve. We also found that the spin-lattice relaxation 
rate in NMR measurement shows a power-law divergence 
1/Ti cx y-i+U/ 2 ^) a t low temperature for any j. 
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APPENDIX: DERIVATION OF CORRELATORS 



We briefly explain the derivation of Eqs. ( |l l| ) and (12). 
As is always the case with the bosonization, we need to 
introduce a short-distance cutoff to obtain finite results. 
The lattice spacing in the original Hamiltonian serves as 
the natural cutoff scale. 

When we use Eqs. (|J) and (^) with the mode expan- 
sions (Q) and (||), we encounter the summation 



oo 

E 1 



1 — COS 



iml 
L + l 



(Al) 



which is formally divergent. We regularize it by inserting 
an exponential factor e - 7r "/( i + 1 ) : 



E ie ^ 



-7m/(L + l) 



1 — cos 



ll 



L + l 



= ln[/(l)], (A2) 



where /(/) = fi(l) defined in Eq. (15). We note that 
Eq. (A.2) is a very good approximation except near the 
points where fi(l) is divergent. Taking derivatives with 
respect to I, we obtain 



-7m/(L + l) 



sin 



71=1 

and 



irnl \ 1 / ttI 

= - cot — r 

L + l) 2 \2{L+1) 



(A3) 



-im/(L+l) 



COS 



iml 



4 sin 



™ 2 



(A4) 



Another point to note is that for Si = ±1 



rl/R 

/ e i27r_Rei0(;) e j27r_R£ 2 0(i')^ ^ / Re i27rR( - ei+e2 ^° d(j>Q 

Jo 

oc (5 ei+£2 , . (A5) 

With the above-mentioned formulas, it is straightfor- 
ward to obtain 



(cos[2tt.R0(O] cos[27ri?0(Z')]) = 



1/(20/(20] 



2[f{l-l')f{l + l')]v 



a i27rRei^(i) -i2nRe!$(l') ie 2 <f>(l' 



ieie 2 e 



S gn(l-l')[f(2l)f(2l')}^ 



-eie 2 e 



hi 



if(i-i>)f(i+ii)nf(2i>)}^> 



l e i2nRe 4>{l) e ieKf>(l)/R e -i27rRe 4>(l') e ie 2 <l>(l')/R 



[f(l-l')f(l + l')}v \f(l + l>) 



1 ( 1 



dl dl' I 2ti \f 2 {l~V) f 2 (l + l')J V2vr 

q \ . M') \ cosiql') g(l + I') - g(l - I') 
sin ■ 



dl 2tt)~"' R I ' 2-kR [/(2/')] 1/2 " 

ie^m/RJe^il'^/R^ _ e / JV_JJ_ \ 

[f(2l)f(2l')]^n \f(l + l') / 



(e :™**-v" ~e 
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TABLE I. The correlation amplitudes; (a) A x = c 2 /2 estimated from the data of (Sf Sf t ); (b) A z = a 2 /2 estimated from 
the data of (Sf). The figures in parentheses indicate the error bar on the last quoted digits. The error bars of A z for A = 
and 0.05 < M < 0.45 are smaller than 10 -5 . The exact values for Mch = given by Eqs. ( po| ) and ( pi] ) are also listed. 
(&)A X 



M 


;h 





0.05 


0.10 


0.15 


0.20 


0.25 


0.30 


0.35 


0.40 


0.45 


A = 


0.0 


0.14709 


0.14626(1) 


0.14364(6) 


0.1390(1) 


0.13262(1) 


0.12410(6) 


0.1132(3) 


0.0993(7) 


0.081(2) 


0.0594(7) 


A = 


0.1 


0.14451 


0.14369(7) 


0.1413(1) 


0.1371(1) 


0.13101(7) 


0.12293(2) 


0.1125(3) 


0.0991(7) 


0.081(2) 


0.0597(7) 


A = 


0.2 


0.14187 


0.1408(4) 


0.1390(2) 


0.1351(2) 


0.1294(1) 


0.12174(5) 


0.1111(7) 


0.0988(6) 


0.081(2) 


0.0600(7) 


A = 


0.3 


0.13921 


0.1384(3) 


0.1366(3) 


0.1330(3) 


0.1278(2) 


0.12053(9) 


0.1111(2) 


0.0985(6) 


0.081(2) 


0.0601(7) 


A = 


0.4 


0.13656 


0.1358(3) 


0.1342(4) 


0.1310(3) 


0.1261(3) 


0.1193(1) 


0.1104(1) 


0.0982(6) 


0.081(2) 


0.0603(7) 


A = 


0.5 


0.13400 


0.1332(4) 


0.1318(5) 


0.1289(4) 


0.1245(3) 


0.1182(2) 


0.10973(9) 


0.0979(5) 


0.081(2) 


0.0605(7) 


A = 


0.6 


0.13164 


0.1310(5) 


0.1294(6) 


0.1268(5) 


0.1229(4) 


0.1170(2) 


0.10905(7) 


0.0976(5) 


0.081(2) 


0.0606(7) 


A = 


0.7 


0.12973 


0.1281(6) 


0.1270(7) 


0.1248(5) 


0.1213(4) 


0.1159(3) 


0.10839(6) 


0.0973(5) 


0.081(2) 


0.0607(7) 


A = 


0.8 


0.12896 


0.1257(8) 


0.1247(8) 


0.1227(6) 


0.1197(5) 


0.1148(3) 


0.10775(7) 


0.0970(5) 


0.081(1) 


0.0609(7) 


A = 


0.9 


0.13214 


0.1233(9) 


0.1223(9) 


0.1207(7) 


0.1182(6) 


0.1137(3) 


0.10714(8) 


0.0967(4) 


0.081(1) 


0.0610(8) 


A = 


1.0 




0.121(1) 


0.120(1) 


0.1188(8) 


0.1177(9) 


0.1127(4) 


0.1065(1) 


0.0958(6) 


0.081(1) 


0.0610(8) 



(b)A z 








0.05 


0.10 


0.15 


0.20 


0.25 


0.30 


0.35 


0.40 


0.45 


A 


= 0.0 


0.05066 


0.05066 


0.05066 


0.05066 


0.05066 


0.05066 


0.05066 


0.05066 


0.05066 


0.05066 


A 


= 0.1 


0.05929 


0.0599(7) 


0.0581(3) 


0.0567(6) 


0.0544(1) 


0.0537(7) 


0.0513(6) 


0.0516(6) 


0.049(1) 


0.0510(8) 


A 


= 0.2 


0.06891 


0.071(1) 


0.0662(6) 


0.063(1) 


0.0580(2) 


0.056(1) 


0.052(1) 


0.052(1) 


0.048(2) 


0.051(1) 


A 


= 0.3 


0.07978 


0.083(2) 


0.0748(7) 


0.069(1) 


0.0614(4) 


0.059(1) 


0.052(1) 


0.053(1) 


0.048(2) 


0.052(2) 


A 


= 0.4 


0.09231 


0.097(3) 


0.0838(6) 


0.075(1) 


0.0645(6) 


0.060(1) 


0.053(2) 


0.053(2) 


0.047(3) 


0.052(2) 


A 


= 0.5 


0.10713 


0.113(5) 


0.093(4) 


0.080(1) 


0.0674(8) 


0.062(2) 


0.053(2) 


0.053(2) 


0.046(3) 


0.052(3) 


A 


= 0.6 


0.12539 


0.132(6) 


0.10263(5) 


0.0854(9) 


0.070(1) 


0.063(2) 


0.054(2) 


0.053(2) 


0.046(3) 


0.052(3) 


A 


= 0.7 


0.14930 


0.153(8) 


0.1121(5) 


0.0903(4) 


0.072(1) 


0.065(2) 


0.054(2) 


0.053(2) 


0.045(4) 


0.052(3) 


A 


= 0.8 


0.18414 


0.176(10) 


0.121(1) 


0.09486(6) 


0.074(2) 


0.066(2) 


0.054(2) 


0.054(2) 


0.045(4) 


0.052(4) 


A 


= 0.9 


0.24844 


0.20(1) 


0.131(2) 


0.0990(7) 


0.076(2) 


0.067(1) 


0.054(2) 


0.054(2) 


0.047(4) 


0.052(4) 


A 


= 1.0 




0.23(1) 


0.139(3) 


0.103(1) 


0.078(2) 


0.067(1) 


0.054(3) 


0.054(2) 


0.044(4) 


0.052(4) 
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FIG. 1. (a) (-l)!'-''!^^) versus r = \l-l'\, (b) |(SfSf,)| 
versus r, (c) (5f ) versus / for A = 0.5 and M ch = 0.05, 0.25, 
and 0.45. The open symbols are the DMRG data and small 
solid symbols are the fits. The numerical errors of the DMRG 
data are smaller than the size of the open symbols. The data 
of (-iy i - l '^{S?S?,) for M ch = 0.05 and 0.45 in figure (a) are 
multiplied by a factor of 2 and 0.5, respectively. 
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FIG. 2. The results of K = 1/(4tt_R 2 ) estimated from the 
fitting of (SfSf,) for < A < 1.0. The dotted curves repre- 
sent the exact values obtained from Bethe ansatz. 
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(a) (-l) 1 '-'''^,,^,,,) versus r = \l - (b) 
\{So,lSo,l')\ versus r , ( c ) (5o,i) versus / for j = 10.0 and the 
Heisenberg case (A = 1.0). The open symbols are the DMRG 
data and small solid symbols are the fits. The numerical er- 
rors of the DMRG data are smaller than the size of the open 
symbols. The data of ( — l)' i_i '' (5™j5™j,) for M = 0.10 and 
0.90 in figure (a) are multiplied by a factor of 2 and 0.5, re- 
spectively. 
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FIG. 4. (a) (-ljI'-'V-,^,,) versus r = \l - (b) 
\{Sg tl Sg ;/}| versus r, (c) (S^,) versus I for j = 1.0 and the 
Heisenberg case (A = 1.0). The open symbols are the DMRG 
data and small solid symbols are the fitting results. The nu- 
merical errors of the DMRG data are smaller than the size 
of the open symbols. The data of (-1) 1 '-'' 1 (S* ,£*,,) for 
M — 0.10 and 0.90 in figure (a) are multiplied by a factor of 
2 and 0.5, respectively. 
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FIG. 6. The correlation function (-l) 1 '"'' 1 {S^^Sl^) in the 
Heisenberg ladder at j — 0.5 and M — 0.2. The dotted curve 
is a guide to the eye. Inset: The M dependence of the IC 
wavenumber q of in the Heisenberg ladder with 

j = 0.5. The dotted line represents the theoretical prediction 
q/TY = M. 
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